Melnikov Analysis for Singularly Perturbed DSII Equation 

Yanguang (Charles) Li 

Abstract. Rigorous Melnikov analysis is accomplished for Davey-Stewartson 
II equation under singular perturbation. Unstable fiber theorem and center- 
stable manifold theorem are established. The fact that the unperturbed ho- 
moclinic orbit, obtained via a Darboux transformation, is a classical solution, 
leads to the conclusion that only local well-posedness is necessary for a Mel- 
nikov measurement. 



1. Introduction 

To build a Melnikov analysis for high dimensional nonlinear wave equations, we 
consider the Davey-Stewartson II equation (DSII) under a singular perturbation 



(1.1) 



iq t = Tq+ [2(\q\ 2 - uj 2 ) + u y ]q + ie[Aq - aq + (3} , 
Au = -4d y \q\ 2 , 



where q is a complex- valued function of the three variables (t,x,y), u is a real- 
valued function of the three variables (t, x, y), the external parameters uj, a, and (3 
are all positive constants, and e > is the perturbation parameter, 

T — 8 — r) A — r) 4- 8 i = \/^-\ 

Periodic boundary condition is imposed, 

q(t,x + 2w/Ki,y) = q(t,x,y) = q{t,x,y + 2tt/k 2 ), 

u(t,x + 2n/ni,y) = u(t,x,y) = u(t,x,y+ 27t/k 2 ), 
where k,\ and «2 are positive constants. Even constraint is also imposed, 
q(t, -x, y) = q(t, x, y) = q(t, x, -y), 

-x, y) = u(t, x, y) = u(t, x, -y). 
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Further constraints are placed upon u>, a, /?, K\, and «2- The first one < aw < (3 
is the condition for the existence of a saddle, and the second one is the condition 
for the existence of only two unstable modes, 



Davey-Stewartson II equation can be regarded as a generalization of the ID cubic 
nonlinear Schrodinger equation (NLS) Q. In fact, it is a nontrivial generalization 
in the sense that the spatial part of the Lax pair of the DSII is a system of two 
first order partial differential equations, for which there is no Floquet discriminant 
to describe the isospectral property, in contrast to the case for NLS. It turns out 
that Melnikov vectors can still be obtained through quadratic products of Bloch 
eigenfunctions, instead of the gradient of the Floquet discriminant as in the NLS 
case. 

At the moment, there is no global well-posedness for DSII in Sobolev spaces. 
In fact, DSII has finite-time blow-up solutions in H S (R 2 ), (0 < s < 1), f§|. Of 
course, DSII has local well-posedness in Sobolev spaces |3| [Q . The Melnikov mea- 
surement is built upon an unperturbed homoclinic orbit of the unperturbed DSII. 
Explicit expression of such a homoclinic orbit can be obtained through Darboux 
transformation The homoclinic orbit is a classical solution. This enables us 
to iterate the local well-posedness result in time, and complete a Melnikov mea- 
surement. Unstable fiber theorem and center-stable manifold theorem are of course 
needed, and established along the same line as in JlJ. Novelties in regularity are 
introduced by the singular perturbation eA which generates the semigroup e €tA . 

The article is organized as follows: section 2 deals with local theory which 
includes unstable fiber theorem and center-stable manifold theorem, and we handle 
global theory in section 3 which includes integrable theory and Melnikov analysis. 



One can view the perturbed DSII as an evolution equation in the q vari- 
able. First, one can define the spatial mean as 



(1.2) 




or 



(1.3) 




2. Local Theory 




q dxdy. 



Then one may introduce the space H s as 



H s = {qeH 



S 



I (?) = o}. 
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2.1. Change of Coordinates. Denote by II the 2D subspace 

(2.2) Tl = {q\d x q = d y q = 0}. 

Dynamics in II is the same as that given in |Jjj|. Denote by S u the circle 

(2.3) S u = {qell\ \q\=u}. 

When acu < (3, there is a saddle Qe near in II, which is located at q = Ie l6 
where 



2.4 I /? 

Its eigenvalues are 

(2.5) ni ti = ±y/i\\Wlpwa6-t 



6(0,f) 



(3 sin t 

~7T 



where / and are given in (2.4). In the entire phase space, Q t is still a saddle. 
Local theory will be built in a tubular neighborhood of S u . Let 

q(t,x,y) = {p(t) + f(t,x,y)}e^ t \ (/) = 0. 

Let 

I=(\q\ 2 )=p 2 + (\f\ 2 ), J = I-u 2 . 



In terms of the new variables (J, 8, /), Equation (2.1) can be rewritten as 
(2.6) j = e 



(2.7) 

(2.8) 
where 



- 2q(J + uj 2 ) + 2/3\/J + uj 2 cos 6 
sin 9 



ell 



2 • 



-2J-e/3 



LJ 



= LJ + VJ - iAf 2 - iMz, 

-*Tf + e(A - a)f - 2na 2 A- 1 T{f + /), 

sinf 



i2JA~ 1 T(f + f) +iefi 



-J, 



-2(V/- V/) +2/3 cos 6* 



-<(/ + /)A- X T(/ + /)) - p-^C/ + /lA-^l/l 2 ) 
-e/Jsinfl' 



A/2 
A/3 



= 2p 



A-^l/l 2 + /A- X T(/ + /) - (fA^Tif + /)> 



= -((/ + /)A- X T(/ + /))/ - 2(|/| 2 )A- 1 T(/ + /) 



/A^TI/I 2 - {/A- 1 !]/] 2 ) 
1 



+2 

— e/3 sin6> 



/9 - 1 ((/ + /)A- 1 T|/| 2 )/ 



VJ + ^ 2 -(|/| 2 ) VJ+^ 
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Since H s (s > 2) is a Banach algebra [Q, we have 

\K J 2 \~0(\\f\\l), i^i-oai/ii;), 
\\Af 2 \\ s ~0(\\f\\ 2 s ), ||A/- 3 ||.~ 0(11/112), (s>2). 



2.2. Unstable Fibers. On II (2.2), the saddle Q e has an unstable and a 
stable curves which lie in an annular neighborhood of S u in II. The width of this 
annular neighborhood is of order 0(y/e). 

Defi nitio n 2.1. Fo r any 8 > 0, we define the annular neighborhood of the 
circle S u (§]§ in II (|J) as 

(2.9) A(S) = {{J,8) | |J| < 8}. 

Unstable fibers with base points in A(8) for some 8 > persist, even under the 
singular perturbation. 

The spectrum of L e consists of only point spectrum. The eigenvalues of L e are: 



(2.10) 



4 



where £ = (6,6), 6 = k j = 0, 1, 2, • ■ ■ , (j = 1, 2), h + k 2 > 0, |£| 2 = £ 2 + £ 2 2 , 
and «i, «2, and w satisfy the constraint (1.2) or (1.3). 

Denote Q ^ by /i^ and /i^ , by /i^. The eigenf unctions corresponding to 
jj^ and fXy are 



e 1 cos K\X 



e ' * cos K 2 y 



2w 

«2 ± l-\/4ti) 2 



2w 



Notice also that the singular perturbation — e|£| 2 breaks the gap between the center 
spectrum and the stable spectrum. Nevertheless, the gap between the unstable 
spectrum and the center spectrum survives. This leads to the following unstable 
fiber theorem. 

Theorem 2.2 (Unstable Fiber Theorem). For any s > 2, there exists a 8 > 
such that for any p G A(8), there is an unstable fiber J 7 ™ which is a 2D surface. Tp 
has the following properties: 

1. Tp is a C 1 smooth surface in || \\ s norm. 

2. is also C 1 smooth in e, a, (3, uj, and p in || \\ s norm, e G [0, eo) for some 
eg > depending on s. 

3. p G , J-p is tangent to span {u^ ,Uy} at p when e = 0. 

4. Tp h as th e exponential decay property: Let S l be the evolution operator of 

\\S l pi - S*p||, < Ce^ +t \\ Pl -p\\ s , Vt < 0, 
where = min{/x+, }. 



5- {3~p } p eA(S) forms an invariant family of unstable fibers, 
S'^c^p, VtG[-T,0], 
and VT > (T can be +oo ), such that S T p G A(8), Vr G [— T, 0] . 



MELNIKOV ANALYSIS FOR SINGULARLY PERTURBED DSII EQUATION 5 

The proof of this theorem follows from the same arguments as in Q. Notice, 
in particular, that C H s for any s > 2. It is this fact that leads to the C 1 
smoothness of in e. Denote by W u (Q e ) the unstable manifold of the saddle 
Q e (|2.4D , which is 3-dimensional. Denote by W^(Q e ) the unstable curve of Q € in 
n ly). Wft{Q e ) = nn W u (Q e ), and Wfi(Q e ) C A(S). W u (Q e ) has the fiber 
representation 

(2.11) W u (Q e )= U 

P eW<f(Q e ) 

Thus W U (Q £ ) C H s for any s > 2. 



2.3. Center-Stable Manifold. Also due to the fact that the gap between 
unstable spectrum and center spectrum survives under the singular perturbation 



(2.1C), a center-stable manifold persists. 



Theorem 2.3 (Center-Stable Manifold Theorem). There exists a C 1 smooth 
codimension 2 locally invariant center-stable manifold W£ s in H n for any n > 2. 

1. At points in the subset W^X^ of W£ s , W^ s is C 1 smooth in e, in H n norm, 
for e E [0, eo) and some €q > 0. 

2. W£ s is C' 1 smooth in (a,/3,w). 



3. The annular neighborhood A(S) in Theorem 2.i is included in W, 



The proof of this theorem follows from the same arguments as in |l| . 
Regularity of W% s in e is crucial in Melnikov analysis. Melnikov integrals are 
the leading order terms in e of the signed distances between W u (Q e ) ( [2.11 ) and 



. The signed distances are set up along an unperturbed homoclinic orbit, and 
the regularity of W^ s in e at e = determines the order of the signed distances in 
e. Due to the singular perturbation, is not C 1 in e at every point rather at 
points in the subset W^ s +A . Here one may be able to replace W^ s +4 by W^+ 2 - But 
we are not interested in sharper results, and the current result is sufficient for our 
purpose. 



2.4. Local Well-Posedness. Following a much easier argument than that in 
c -g- @ @li one can prove the following local well-posedness theorem. 

Theorem 2.4. For any q € H n (n > 2), there exists t = r(||g ||n) > 0, 
such that the perturbed DSII (2.1) has a unique solution q(t) — S (qo;e,a,f3,u>) S 
C°([0, t], H n ), q(0) = qo, where 5* denotes the evolution operator. S 1 ^-; e, a, j3, w) : 
H n i ^ H n is C 1 in q and (a,/3,u). 5*(-; e, a, (3, w) : H n+i ^ H n is C 1 in t and 
e, e £ [0,e ), e > 0. 

Here "C 1 in q and (a,/3,w)" can be replaced by "C°° in q and {a,j3,u)". 
H n+4 can be replaced by H n+2 . But we are not interested in sharper results. 



3. Global Theory 

Global theory is referred to a theory global in phase space, which includes 
integrable theory and Melnikov analysis. Integrable theory provides two ingredients 
for a Melnikov analysis: (1). An explicit expression of the unperturbed homoclinic 
orbit, (2). Melnikov vectors with explicit expressions. 
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3.1. Integrable Theory. Calculations in this subsection are essentially the 
same with those in |^| . The minor differences are introduced by the spatial periods 
2tt/ki and 2ir / k 2 in contrast to 2tt and 27r in 0]. Proofs of theorems and lemmas 
can be found in |^| . 



The DSII [e = in (1.1)] is an integrable system with the Lax pair 

(3.1) Lip = Xtp , 

(3.2) d t ip 



A,p, 



where ip = (ipi,ip2) , and 



D- 



A 



-d 2 x qd x 
qd x dl 



where 

(3.3) D+ = ad y + d x , 

t*i and r 2 have the expressions, 



D~ 



ri (D+q) 
-{D-q) r 2 



ady - d x , 



(3.4) 



1 



n = -[-2(| 9 | 2 - uj 2 ) - u y + iu] , r 2 = -[2(\q\ 2 - u 2 ) + u y + li 



where u is also a real-valued function satisfying 

Ait = Aiadxdy^ 1 ^ . 
Notice that DSII is invariant under the transformation a: 
(3.5) a(q, q, n, r 2 ; a) = (q, q, -r 2) -n; -a) 



Applying the transformation <j (|3.5| ) to the Lax pair (3^, ^j), we have a congruent 
Lax pair for which the compatibility condition gives the same DSII. The congruent 
Lax pair is given as: 



(3.6) 
(3.7) 

where tp — (ipi,ip2), and 



L = 



Lip — Xip , 



<1 

D- 




A 



-d 2 x qd x 
qd x dl 



(D+q) 



-{D-q) 
-ri 



The Backlund-Darboux transformation can be formulated as follows. Let (q, u) 
be a solution to the DSII, and let A* be any value of A. Let ip — (ipi,4>2) T be a 
solution to the Lax pair ( |3.l| , 3.2 ) at (q, q, ri, r 2 ; A*). Define the matrix operator: 

A + a b 
c A + d 



r = 
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where A = ad y — A, and a, b, c, d are functions defined as: 



[lp2 A 2 lp2 
[i>2 Al -01 



d = 



— [-01 Al 1p2 

— \jp2 Al 1p2 



■ 4>i Ai ipi] 

4>2 A 2 -0l] 

■ ipi A 2 -0il 



in which Ai = ad v — A*, A 2 = a<9j, + A*, and 

A = -[|^| 2 + h/> 2 | 2 
Define a transformation as follows: 

(q,n,r 2 ) - 



(3.8) 



(Q, Ri,Ri) , 



Q = q-2b, 
i?i = n + 2(D + a) , 
i? 2 = r 2 -2(D-d), 



where <p is any solution to the Lax pair (3.1, p.2| ) at (q, q, n, r 2 ; A), D + and I? are 
defined in (|3.3|), we have the following theorem [H]. 



Theorem 3.1. The transformation (\3.8\) is a Backlund-Darboux transforma- 
tion. That is, the function Q defined through the transformation (3.8) is also a so- 
lution to the DSII. The function $ defined through the transformation fc>.8\ ) solves 
the Lax pair §T% filfy at (Q, Q, R 1 ,R 2 ;X). 

Consider the spatially independent solution, 

(3.9) q c = r/cxp{— 2i[q 2 — uj 2 ]t + ij} , 



where rj satisfies the constraint (1.2) and (1.3) with u> replaced by r\. The dispersion 
relation for the linearized DSII at q c is 



n = ± 



V - (£? + &) , for 5q ~ q c exp{i(^x + &») + fit} , 



where £i = fciKi, £ 2 = fc 2 K 2 , and k\ and fc 2 are integers. There are only two 
unstable modes (ki,0) and (0, k 2 ) under even con straint. 



The Bloch eigenfunction of the Lax pair (|3.l[ ) and (3.2) is given as, 



(3.10) 



xjj = c(t) 



x 



exp {i(^x + £ 2 y)} , 
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where 

c(t) = c cxp {[2£i(ia£> - A) + ir 2 ] t} , 
T2-r\ = 2(\q c \ 2 - u 2 ) , 

X = - A) - , 

For the iteration of the Backlund-Darboux transformations, one needs two sets of 

2 — l K i (f° r a fixed 



eigenf unctions. First, we choose £1 = ±|ki, £2 = 0, A = \/ .7 
branch) , 

-9c 



(3.11) 
where 



<\\p <j ±i-KlX 



= CQexp{[=FKiA + ir 2 ]t} 



X =-A T«2 Kl 



1 



-Ki ± «Aq . 



We apply the Backlund-Darboux transformations with ^ = V J+ + V* i which gener- 
ates the unstable foliation associated with the (/ti, 0) linearly unstable mode. Then, 

we choose £2 = ±\k 2 , A = 0, = wry 2 — -jK 2 (for a fixed branch), 



(3.12) 
where 



c± 



-9c 
X± 



cxp \ i(t$x± -K 2 y) } , 



c± = c± exp { [±ia«2^i + ^2] *} 



X ± = ±ia-K 2 - itf = ±ve™ 2 



i.e. rye 



±«2 



1 



ia-n 2 ± z£j . 



We start from these eigenfunctions <j>± to generate T<J)± through Backlund-Darboux 
transformations, and then iterate the Backlund-Darboux transformations with T<p + + 
to generate the unstable foliation associated with all the linearly unstable 
modes (k\,0) and (0,n 2 ). It turns out that the following representations are con- 
venient, 



(3.13) 

(3.14) 
where 



\J c o c o 



,ir 2 t 



"1 



«r = 



± = rye^* 5 



and 



w, = 



w 2 = ±rye 2 



T = 2n 1 \ t-p, x = -k 1 x + -, z 
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/„0 _ e p+i# 



t = 2iaK 2 £°t + p , y=-K 2 y+-, z = y - d 



2t ( vi 

V 2 



The following representations are also very useful, 

(3.15) ^ = ^ + + ^- = 2^4^^( 

(3.16) <(> = 4> + + r ^^4^- e<x+ir2t {Z 2 ) ^ 

where 

T T T T 

V\ = —q c [cosh — cos x — i sinh — sin x] , v 2 = r\ [cosh — cos z — i sinh — sin z ] , 
io i = — g c [cosh — cosy + i sinh — sin y] , W2 = Tj [sinh — cos z + i cosh — sin z] . 



Applying the Backlund-Darboux transformations ( |3.8[ ) with -0 given in (3.15), we 
have the representations, 



(3.17) 



b = -q c b = - 



Xq sech t sin(5 + z) sirt(x — z) 

1 + sech r cos(a; + z) cos(x — z) 
Ao<7< 



cos(i — z) — i tanh r sin(i — z ) 



1 + sech t cos(x + z) cos(a; — 2) 



(3.18) + sech r cos(i + z) 

(3.19) c = b, d = —a=—a. 

The evenness of b in x is enforced by the requirement that i? — #i = ±-|, and 



(3.20) 



=FAo sech rcos^i sin(Kix) 
x 1 =p sech r sin i?i cos (kix) 



V 



— sin i?i — 4 tanh r cos $1 



1 =F sech t sin #1 cos(«ia;) 



(3.21) ± sech tcos(kix) 

(3.22) c=b, d = -a=-a. 

Notice also that a ± is an odd function in x. Under the above Backlund-Darboux 



transformations, the eigenfunctions <j>± ( 3.12j ) and <fi are transformed into 



(3.23) 
where 



ip = r</>± , ip = r<j) = r<f>+ + !>_ , 



r = 



A + a b 
b A - a 
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and A = ad v — A with A evaluated at 0. Then 



where 



-QcW? 



Wf 



wf 



[±i-an 2 +a± rjbe^^i^ , 
± e ^[±i-aK 2 -a± j^e ± ** a ]e ± * ±il > 



-q c W 1 
r]W 2 



where 



f 1 
W\ = cosh — [a cos y — — olk% sin y + ir/b sin z] 

t 1 . 

+ sinh — [—iceK2 cos y + ia sin y + ijb cos f ] , 

f . . A 1 . 
= cosh — [— zasinz + —iaK2 cosz + rjb cosy] 

f 1 . A = 

+ sinh — [— — O.K2 sin z — a cos z + irjb sin y] . 



We generate the coefficients in the Backlund-Darboux transformations ( |3.8| ) with 
tp (the iteration of the Backlund-Darboux transformations), 



(3.24) a (7) 

(3.25) 
(3.26) 



9c 



|Wi| 2 + |W 2 | : 



n -1 



|Wa| 2 + IW2I 2 



3« = bW , = - a M , 
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where 



W 2 (ad y W 2 ) + W 1 (ad y W 1 ) 



—aK2< coshf 



2 2 2|t|2 

4*2 -a -rj \b\ 



|Wil 2 + |VF 2 I 2 



an 2 a + ian(b + b) cos d 2 
cos(y + 5) sini9 2 + sinhf 



arj(b — b) sin $2 



= cosh t 
1 



fl2 + ^ K 2 + V 2 \b\ 2 + iaK2i]-(b + b) cos #2 



2 „2 „2lr|2 



sin(y + z) sin t?2 + sinh f 



C^K 2 n- ~ ^) sm ^2 



W 2 (a^Wi)-Wi(^W 2 ) 



aK 2 ?76 + i(-a 2 + -k| + t? 2 ^ 2 ) cos tf 2 



sinh f 



2 2 1 2I2 
a - T K 2 + ?7 



sini?c 



The new solution to the DSII is given by 

(3.27) Q = q c - 2b- 2b {1) . 

The evenness in y is enforced by the requirement that {} — $2 = i?. In fact, 

we have 

Lemma 3.2. Choosing the Backlund parameters •& and-d as follows: i? = $!±-| ; 
and = tf 2 ± § , 

(3.28) &(-!)= 6(1) , &W(-a;,i/) = 6«(a; ) i/)=6C J )(a: > - 1 /), 

and Q = q c — 2b — 2b^ is even in both x and y. 

The asymptotic behavior of Q can be computed directly. In fact, we have the 
asymptotic phase shift lemma. 

Lemma 3.3 (Asymptotic Phase Shift Lemma). For A > 0, £° > 0, and a — 

—i; as t — > ±00, 

(3.29) Q = q c -2b- 26« q.e^e^ 2 ^ 1 '^ . 

In comparison, the asymptotic phase shift of the first application of the Bdcklund- 
Darboux transformations is given by 

q c -2b^ q c e^ . 



Next we generate the Melnikov vectors. Starting from ip^ and <j>± given in (3.11) 
and (3.12), we generate the following eigenfunctions corresponding to the solution 
Q given in (3.27) through the iterated Backlund-Darboux transformations, 



(3.30) 
(3.31) 



V* = r^IV*, at X = \ = \jn 2 ~ 1 K 2 1 

$± = r (/) r</>±, atA = o, 
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where 



r = 



A + a b 



b A - a 



A + a« 



where A = ad v — A for general A. 

Lemma 3.4. The eigenf unctions and $>± defined in ( 3. Si ) and ( 3.3\ ) have 
the representations, 



±i\ K 1 r ) - 1 JcJc^e ir * t [\v 1 \ 2 + \v 2 \ 2 }- 
-q c (A - aW)W + ^ (J) W 



(3.32) 



- T]bWv2- (A + a^)n 
(3.33) *± = ±i~aK 2 Jc^e i £ x ^ r2t [\Wi\ 2 + \W 2 \ 2 }- 1 



-Qc^i 

V^2 



where — —qjM\ and 

Ex = 2W[{WfWi) + W 2 T {W^W^) + Wf(W 1 ~W 2 + ) , 

s 2 = 2W 2 {w^w^) + wf{w}w^ ) + wf {w^wf). 

If we take r 2 to be real [in the Melnikov vectors, r 2 appears in the form r 2 — r\ 
2{\q c \ 2 -u 2 )], then 



(3.34) 



-> , $± -> , ast^ ±oo . 



Next we generate eigenfunctions solving the corresponding congruent Lax pair 
( 3.6 , 3.7 ) with the potential Q, through th e ite rated Backlund-Darboux transfor- 
mations and the symmetry transformation (3.5) J5). 

Lemma 3.5. Under the replacement 

a — ► —a (then d 2 — > 7r — "d 2 ), T\ — > —r 2 , 

(3.35) r 2 — >- ri , d — >d + n-2d 2 , p — ► -p, 

the potentials are transformed as follows, 

Ri — ► — R 2 , 
R 2 — ► —R\ . 



The eigenfunctions , 3> ± and <&± given in ( |02j ) and ( |3.33| ) depend on the vari- 
ables in the replacement (3.35): 



*± = 
$+ = 



* ± (a,r 1 ,r 2 ,i?, i o). 
$ ± (a,r 1 ,r 2 ,d,p) . 
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Under replacement (3.35|), , J ± and <i>± are transformed into 



1.3 



(3.36) 
(3.37) 



$ ± = * ± (-a,-r 2 ,-ri,^ + 7r-2i?2 ) -p) 
$± = $±(-Q!,-r2,-ri,i9 + 7r-2i?2,-p). 



Lemma 3.6. and$± solve the congruent Lax pair (3J^ 3. 7) at (Q, Q, R±, i? 2 ; Xq) 
and (Q, Q, Ri,R2] 0), respectively. 

Notice that as a function of 77, £® has two (plus and minus) branches. In 
order to construct Melnikov vectors, we need to study the effect of the replacement 

Lemma 3.7. Under the replacement 
(3.38) £ — >-£ (then #2 — > -0 2 ), 3 — >tf + 7r-2tf 2 , p — > -p, 
the potentials are invariant, 

Q — ► Q , Ri ► Ri , i?2 ► R2 ■ 



The eigenfunction <&± given in ( |3.33 ) depends on the variables in the replace- 
ment (|3~38l): 



Under the replacement ( |3.38] ), $± is transformed into 
(3.39) *± = *±(-€?,^ + 7r-2i? 2> -p) 



Lemma 3.8. $± sofces tfte Lax pair (3.1,3.i) at (Q, Q, Ri, R 2 ; 0) 



In the construction of the Melnikov vectors, we need to replace $± by $± to 
guarantee the periodicity in x of period 

The Melnikov vectors for the Davey-Stewartson II equations are given by, 



(3.40) 



*2~*2 
* + $ + 



5 



4-2^ 



(3.41) 



U 4 



$(2)^(2) 



$(2) $ (2) 



where "-" denotes complex conjugate, and S = 



1 

1 



In fact, the even parts 



of U + and U+ are the Melnikov vectors in our phase space. Nevertheless, the 
Melnikov integral formulas end up the same, as shown in J5|. For simplicity, we 
just use U + and U+. 
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3.2. Melnikov Analysis. The main difficulty in a rigorous Melnikov mea- 
surement is due to the lack of global well-posedness. The main idea in resolving 
this difficulty is to iterate the small time interval in local well-posedness by virtue 
of the fact that the unperturbed homoclinic orbit is a classical solution. 

Let p be any point on W^(Q e ), the unstable curve of Q e in IT. By the Unstable 
Fiber Theorem fT2[ is C 1 in e for e G [0, eo)> £o > 0; thus, there are two points 
q e (0) and <7o(0) on the unstable fibers and T£ | £= o, such that 

WqM-Uomin^C^e, (n>2). 
The key point here is that C H s for any fixed s > 2. The expression of the 



unperturbed homoclinic orbit qo{t) has been given in (3.27) which represents a 
classical solution to the DSII. Let 

D* s = sup { \\q (t)\\s } , 0>2). 

£6 ( — oo.+oo) 



By the Local Well-Posedness Theorem |2.4[ there exists r — r(Z?* ) > 0, such that 

he{t) - qo(t)\\n < C^e , t€[0,T] , 
where = Cn 2 ^(D* +4 ). There is an integer N > such that 

gb(JVr) G W n cs |e=o , 



where W£ s is given by the Center-Stable Manifold Theorem |2.3| . Iterating the Local 
Well-Posedness Theorem AT times, one gets 

\\qe(t) - qo(t)\\n <C^e , tE[0,NT] , 

where = Cn^ (D* , 4 ). Our goal is to determine when q e (Nr) G W^f through 



Melnikov measurement. The two Melnikov vectors U + and U+ ( |3.40| )-(3.41) are 
transversal to W£ s - There is a unique point q e (Nr) G W£ s such that 

q e (Nr) - q e (N T ) G span {U + \U+} 5 



thus, q e {Nr) G W^ 4 . By the Center-Stable Manifold Theorem |2. 3, 

||,7 e (7Vr)-go(^r)|| n <C( 4 )e, 

where d 4) = ck 4) (D* n+4 ). Thus 

||g e (iVT)-g e (iVr)|| n <C n e , 

where C„ = C n (£)* +4 ). To determine when q e (Nr) = q e (Nr), one can define the 
signed distances 

d x = (U+, q e (Nr) - £{Nt)) , d 2 = (U+, q^Nr) - £(Nr)) , 

where q = (q, q) T , and 

(A,B) = / / {Aj\Bi dxdy . 

Jo Jo 

The rest of the derivation for Melnikov integrals is completely standard. For details, 



see e.g. fl9| [10 



d k =eM k + o(e), k = l,2, 

where 



/OO P OO 

(U + ,G) dt , M 2 = / (U+,G) dt , 
-OO </ —OO 
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where G = (/, f) T , f = AQ - aQ + [3. That is, 

/OO p2-K I «2 /»2tT /Ki 
/ / Re {(*l$t)f + (*t*t)f} dxdydt , 
-oo JO JO 

/OO p27T I «2 /•2tt 

/ / Re {($i 2) $f)/ + (3? drdj/dt , 
-oo Jo JO 

where 77 = w, and we divide by the constant zAqKi a/cJcq e 17 / 2 , and $ + by 



jiaK2i]<\J c° + c®_e %1 1 2 . It has been verified numerically that multiplication of \& + and 
<&+ by a complex constant leads to equivalent results. It turns out that 

Mj = Mf ] + aMf ] + (3 cos 7 m] 3) + (3 sin 7 Afj 4) , (j = 1, 2) , 

where Afj = Afj l) (w,Ap), (j = 1,2; 1 < / < 4), Ap = p + iaK^K^Xo l p, 
t = iaK2S,iK^ 1 Xq 1 t + Ap. 

Mj = (j = 1,2) imply that 

a = A/5,7) = jMf [cos 7M 2 (3) + sm 7 M 2 (4 } ] 
-Af 2 (1) [cos 7 Af 1 (3) + sin7M 1 (4) ] 
x { M { 2 2) [cos 7 Af 1 (3) + sin 7 M 1 (4) ] 



(3.42) -Mf> [cos 7 M 2 (3) + sin 7 M^ 4) ] 
[3 = /3(w, Ap, 7) = [M^M™ - Af 2 (1) Aff] 

x |ikf 1 (2) [cos 7 Af 2 (3) + sin 7 A/4 4) 

(3.43) -M 2 (2) [cos 7 Af 1 (3) + sin 7 Af 1 (4) ] 



-1 



Theorem 3.9. There exists €q > 0, swc/i that for any e G (0, eo), iftere exists 



a domain T> e C M x R + x R + where uj satisfies the constraint (l.k) or (l.S), and 



auj < [3. For any (a, (3, of) € T> t , there exists another orbit in W u (Q e ) n W£ s other 



than the unstable curve W^(Q e ) of Q t in II, for the perturbed DSII (1.1). 



Proof. The zeros of Mj (j = 1, 2) are given by ( |3.42| ) and ( |3.43|) . We need a > 



and /3 > which define a region in the external parameter space, parametrized by 
Ap and 7. Then the theorem follows from the implicit function theorem. Q.E.D. 
For example, when K\ = 1 and k 2 = y2) 

a(^ + 0.11, 1.1, |) = 5.645 , /3(^ + 0.11, 1.1, |) = 11.336 . 
4. Appendix 

The main obstacle toward proving the existence of a homoclinic orbit for the 
perturbed DSII ((□]) comes from a technical difficulty in the normal form transform 
Q|. In this appendix, we will present the difficulty. 
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4.1. The Technical Difficulty in the Normal Form Transform. To lo- 
cate a homoclinic orbit to Q e , we need to estimate the size of the local stable 
manifold of Q e . The size of the variable J is of order 0(^/e). The size of the variable 
9 is of order 0(1). To be able to track a homoclinic orbit, we need the size of the 
variable / to be of order 0(e tl ), \i < 1. Such an estimate can be achieved, if the 
quadratic term A/jj in (2.8) can be removed through a normal form transformation. 
In fact, it is enough to remove its leading order part 



A/2 = 2u> 



A^Tljf + /A- 1 T(/ + /) - (/A- X T(/ + /)) 



That is, our goal is to find a normal form transform g = f + K (/, /) where K is a 
bilinear form, that transforms the equation 

ft = LJ - iJ\f 2 , 

into an equation with a cubic nonlinearity 

g t =L e g + 0(\\gf s ), (s > 2), 

where L t is given in (|2.8|). In terms of Fourier transforms, 



fe^O fc^O 

where k = (fci, fe) G Z 2 , £ = {n\x, K2J/). The terms in A/2 can be written as 
A-^l/l 2 = i E a(fc + £) 



f{k)f{-l) + f{t)f{-k) 
fA-'Tf - (/A _1 T/) = i E f a (*) + *m(k)f(i)e l(k+eH , 

k+l^O 

\ E [aW/>)/H)+a(fc)/>)7Pfc) 
L 



where 



k+£=£0 



i(k+t)-t 



a{k) 



&2^2 



L.2 yJl j_ L,2 ^2 

1 1 2 2 



We will search for a normal form transform of the general form, 

g = f + K(fJ), 



where 



K(f,.f) = J2 [Ki(k,i)f(k)f(£) + k 2 (k,l).f(k)f(-e) 



k+e^o 



-K 2 (£,k)f(-k)f(£) + K 3 (k,£)f(-k)f(-£) 



i(k+e)a 



where Kj(k,£), (j = 1,2,3) are the unknown coefficients to be determined, and 
Kj(k,£) — Kj(£,k), (j = 1,3). To eliminate the quadratic terms, we first need to 
set 

iL e K(f, /) - iK(LJ, /) - iK(f, LJ) = A? 2 , 
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which takes the explicit form: 



(o-! + ia)Ki(k, £) + B(£)K 2 (k, £) + B(k)K 2 (£, k) 



(4.1) 

(4.2) 

(4.3) 
(4.4) 



-B(k + £)K 3 (k, £) = —[B(k) + B{£% 

ZtU 



-B(l)Ki(k, I) + (era + icr)K 2 (k, £) + B(k + £)K 2 (£, k) 
+B(k)K 3 (k, £) = ±-[B(k + £) + B(£)}, 



-B(k)K!(k,£) + B{k + £)K 2 (k,£) + (<j 3 + ia)K 2 (£,k) 
+B(£)K 3 (k, £) = ^[B{k + £) + B(k)], 



B{k + £)Ki(k, £) ~ B(k)K 2 (k, £) - B(£)K 2 (£, k) 
+(a 4 + ia)K 3 (k,£) = 0, 



where B(k) — 2ui 2 a(k), and 

(i = 6 a - 2{ki£in\ + k 2 £ 2 K\) 

o-i = 2{k 2 £ 2 n\ - kihul) + B(k + £) - B(k) - B{£) , 
<T2 = 2[(fc 2 + £ 2 )£ 2 k 2 2 - (fci + li)£n4] + B{k +£)- B{k) + B{£) , 
<j 3 = 2[(fc 2 + £ 2 )k 2 n 2 2 - (ki + lx)k X Kl\ + B[k + £) + B(k) ~ B{£) , 
0-4 = 2[(k% + k 2 £ 2 + £\)k\ - (k\ + Mi + £\)n\] + B(k + £) + B{k) + B(£) . 
Since these coefficients are even in (k,£), we will search for even solutions, i.e. 
K j (-k,-£)=K j (k,£), j = 1,2,3. 

The technical difficulty in the normal form transform comes from not being 
able to answer the following two questions in solving the linear system (4.1)-(4.4): 

1. Is it true that for all k,£ € Z 2 /{0}, there is a solution ? 

2. What is the asymptotic behavior of the solution as k and/or £ — > oo ? In 
particular, is the asymptotic behavior like k~ m and/or £~ m (m > 0) ? 

4.2. A Formal Calculation. Formally conducting the calculation for the 
second measurement to locate a homoclinic orbit Q, one gets the formulas 



Mj = (j = 1,2), /?cos 7 =- 



2 sin 



A7 



where A7 = — 4(i9i — $2)- Thus we have a — \ j\ 



-cA 7 [2sin^]' 1 (M 1 (3) M 2 (4) - M^M^) 



(3 = Ap) 



(awA 7 ) 2 [2 sin 



A 7l-J 



+ (M 1 (4) )- 2 [^i (1) + a(Mi (2) - M 1 (3) wA 7 (2 



am—) 



l\l2 



-,1/2 
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For example, when K\ = 1 and k 2 = V%, 

/2 

x{\ + 0.11, 1-1) = 0.4326 . 
References 

[1] Y. Li, Persistent homoclinic orbits for nonlinear Schrodinger equation under singular pertur- 
bation, Submitted (2001). 

[2] T. Ozawa, Exact blow-up solutions to the Cauchy problem for the Davey-Stewartson systems, 
Proc. R. Soc. Lond. 436 (1992), 345. 

[3] J. M. Ghidaglia and J. C. Saut, On the initial value problem for the Davey-Stewartson 
systems, Nonlinearity 3 (1990), 475. 

[4] L. Y. Sung, An inverse scattering transform for the Davey-Stewartson II equations, part I, 
II, III, J. Math. Anal. Appl. 183 (1994), 121, 289, 477. 

[5] Y. Li, Backlund-Darboux transformations and Melnikov analysis for Davey-Stewartson II 
equations, Journal of Nonlinear Sciences 10, no.l (2000), 103. 

[6] R. Adams, Sobolev Space, Academic Press, New York, 1975. 

[7] T. Kato, Nonstationary flows of viscous and ideal fluids in R 3 , Journal of Functional Analysis 
9 (1972), 296. 

[8] T. Kato, Quasi-linear equations of evolution, with applications to partial differential equa- 
tions, Lecture Notes in Math., Springer 448 (1975), 25. 

[9] Y. Li and D. W. McLaughlin, Homoclinic orbits and chaos in discretized perturbed NLS 
system, part I. homoclinic orbits, Journal of Nonlinear Sciences 7 (1997), 211. 
[10] Y. Li et al., Persistent homoclinic orbits for a perturbed nonlinear Schrodinger equation, 
Comm. Pure Appl. Math. XLIX (1996), 1175. 

Department of Mathematics, University of Missouri, Columbia, MO 65211 
E-mail address: cliamath.missouri.edu 



